nature astronomy 


Article 


Q 


https://doi.org/10.1038/s41550-024-02249-0 


Identification of the weak-to-strong 
transition in Alfvénic turbulence from 


space plasma 


Received: 17 January 2023 


Accepted: 11 March 2024 


Siqi Zhao ©”, Huirong Yan®*® 


, Terry Z. Liu®?, Ka Ho Yuen ®° & Huizi Wang ®* 


Published online: 12 April 2024 


® Check for updates 


Plasma turbulence is a ubiquitous dynamical process that transfers energy 
across many spatial and temporal scales in astrophysical and space plasma 
systems. Although the theory of anisotropic magnetohydrodynamic 


(MHD) turbulence has successfully described natural phenomena, 

its core prediction of an Alfvénic transition from weak to strong MHD 
turbulence when energy cascades from large to small scales has not been 
observationally confirmed. Here we report evidence for the Alfvénic 
weak-to-strong transition in small-amplitude, turbulent MHD fluctuations 
in Earth’s magnetosheath using data from four Cluster spacecraft. Our 
observations demonstrate the universal existence of strong turbulence 
accompanied by weak turbulent fluctuations on large scales. Moreover, we 
find that the nonlinear interactions of MHD turbulence are crucial to the 
energy cascade, as they broaden the cascade directions and fluctuating 
frequencies. The observed connection between weak and strong MHD 
turbulence systems may be present in many astrophysical environments, 
suchas star formation, energetic particle transport, turbulent dynamos, 
and solar corona or solar wind heating. 


The theory of anisotropic magnetohydrodynamic (MHD) turbulence 
has been widely accepted and adopted for plasma systems ranging from 
clusters of galaxies, the interstellar medium and accretion disks to the 
heliosphere’*. One of the most crucial predictions of the theory is an 
Alfvénic transition from weak to strong MHD turbulence when energy 
cascades from large to small scales**. The self-organized process from 
weak to strong MHD turbulence is the cornerstone of understanding 
the energy cascade in the complete picture of MHD turbulence. 

The critical balance model is attractive for describing the phys- 
ical behaviour of incompressible MHD turbulence*. When T4 < Tni 
(referred to as weak MHD turbulence), weak interactions among 
the counter-propagating wave packets transfer energy to higher k, 
whereas no energy cascades to higher k,, where T, = 1/(k,V,) is the 
linear Alfvén wave time, T,, = 1/(k, 5V,) is the nonlinear time, V, is the 


Alfvén speed, 8V; is the fluctuating velocity perpendicular to the 
background magnetic field (Bo), and k, and k, are wavenumbers per- 
pendicular and parallel to B, (ref. 6). As turbulence cascades to smaller 
scales, nonlinearity strengthens until it reaches the critical balance 
(Ta = Ta) at the transition scale (A;,). On scales smaller than Acs, Alfvén 
wave packets are statistically destroyed in one T,. In addition to the 
first-order interactions of counter-propagating waves, all higher 
orders of interactions can contribute to create strong MHD turbu- 
lence*”. In compressible MHD, small-amplitude fluctuations can be 
decomposed into three eigenmodes (namely, Alfvén, fast and slow 
modes) ina homogeneous plasma® ”. Alfvén modes decoupled from 
compressible MHD turbulence, which are linearly independent of fast 
and slow modes, show similar properties to those in incompressible 
MHD turbulence, for example, the Kolmogorov spectrum and the 
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scale-dependent anisotropy”. Additionally, numerical simulations 
have confirmed that the Alfvénic weak-to-strong transition occurs in 
bothincompressible MHD turbulence and Alfvén modes decomposed 
from compressible MHD turbulence”, 

However, this transition has not been confirmed from observa- 
tions. Inthis study, we present evidence for the Alfvénic weak-to-strong 
transition and estimate the transition scale Acs in Earth’s magne- 
tosheath using data from the four Cluster spacecraft’. Earth’s mag- 
netosheath offers a representative environment for studying plasma 
turbulence, given that most astrophysical and space plasmas with 
finite plasma £ are compressible, where £ is the ratio of the plasma to 
magnetic pressure. 


Results 

Here we present an overview of fluctuations observed by Cluster-1 
in geocentric-solar-ecliptic (GSE) coordinates during 23:00-10:00 
Universal Time (UT) on 2-3 December in Fig. 1. During this period, 
the four Cluster spacecraft flew in a tetrahedral formation with rela- 
tive separation d.. = 200 km (around three proton inertial lengths 
d, = 74 km) on the flank of Earth’s magnetosheath around [1.2, 18.2, 
-5.7] R,; (Earth radius). We chose this time interval to study the Alfvénic 
weak-to-strong transition because the fluctuations satisfy the follow- 
ing criteria. First, the background magnetic field (B) measured by the 
fluxgate magnetometer” and the proton bulk velocity (V,) measured 
by the Cluster lon Spectrometer (CIS)'* were relatively stable, as shown 
in Fig. 1a,b. Second, we cross-verified the reliability of the plasma data 
based on the consistency between the proton density (N,) measured 
by CIS and the electron density (N.) measured by the Waves of High 
frequency and Sounder for Probing of Electron density by Relaxation 
(WHISPER)”, as shown in Fig. 1c. 

We set a moving time window with a length of 5 h and a moving 
step of 5 min. The selection of a5 h length ensured that we obtained 
measurements at low frequencies (large scales) while the mean mag- 
netic field (B) within the moving time window was approaching the 
local mean field at the selected largest spatial scale. The uniform- 
ity of By) was independent of the transformation between real and 
wavevector space; however, it differed from the theoretically expected 
local mean field at each scale. To assess such differences, Supple- 
mentary Fig. 1 shows that the local mean field at different scales is 
closely aligned with B, most of the time, suggesting that using B, to 
approximate the local mean field is acceptable. To further address 
this limitation of the mode decomposition method, which relies ona 
perturbative treatment of fluctuations in the presence of a uniform 
background magnetic field®, we provide results obtained using vari- 
ous time window lengths, all of which show similar conclusions (Sup- 
plementary Fig. 2). 

Third, Fig. 1d shows the spectral slopes of the trace magnetic field 
and proton velocity power calculated by a fast Fourier transform with 
three-point centred smoothing in each time window. These spectral 
slopes at the spacecraft-frame frequency f,, ~ [0.001 Hz, 0.1f,,] are close 
to -5/3 or -3/2 (the proton gyro-frequency f,; = 0.24 Hz), suggesting 
that the turbulent fluctuations were in a fully developed state. The 
remaining magnetosheath fluctuations with spectra close to fz! are 
typically populated by uncorrelated fluctuations” and are beyond 
the scope of the present paper. Figure le shows that the average proton 
plasma £, was around 1.4. Finally, Fig. 1f shows that the turbulent Alfvén 
number M4 turb = V/V, = 5B/(2B,) = 0.33, suggesting that the fluctua- 
tions include substantial Alfvénic components and satisfy the 
small-amplitude fluctuation assumption (the nonlinear terms ôV? and 
5B’ are weaker than the linear terms V,5V, and BSB). Nevertheless, the 

|6By (foc)? x 0.34, 
|6By (fsc)|2+16B1 (fee) |? 
indicating that the fluctuations are a mixture of Alfvén and compress- 
ible magnetosonic modes (fast and slow)”, where 6B, and 6B, are the 
fluctuating magnetic fields parallel and perpendicular to B,. 
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Fig. 1| Overview of fluctuations measured by Cluster-1 in Earth’s 
magnetosheath on 2-3 December 2003. a-f, Data are displayed in GSE 
coordinates. a, Magnetic field components (Bx, By and B7). b, Proton bulk velocity 
(Vx Vyand V,). c, Proton and electron densities. d, Spectral slopes (a) of magnetic 
field and proton velocity fluctuations between 0.001 Hz and 0.1f,;. The two 
horizontal lines represent a = -5/3 and -3/2. e, Proton plasma £,. f, Turbulent 
Alfvén Mach number (M4 tur» = 5V,/Va) and half of the relative amplitudes of the 
magnetic field (6B/(2B,)), where 5V, and 6B are the root mean square of the 
proton velocity and magnetic field fluctuations, respectively. The fluctuations 
analysed in detail were observed during 23:00-10:00 UT on 2-3 December and 
are marked between the two vertical dashed lines. 


Due to the homogeneous and stationary state of the turbulence 
(Supplementary Fig. 3), we can utilize frequency-wavenumber distri- 
butions of Alfvénic power, that is magnetic power Pg, (k,,k).fgc) and 
proton velocity power Py, (k,,kj.fsc), to investigate the structure of 
the turbulence. Alfvénic fluctuations were extracted in each time 
window based on their incompressibility and fluctuating directions 
perpendicular to B (Methods). To distinguish between spatial evolu- 
tion and temporal evolution without relying on any spatio-temporal 
hypothesis, we determined the wavevectors by combining the singular 
value decomposition method” (to obtain ksyp) and a multispacecraft 
timing analysis” (to obtain k,). Note that k, was not completely 
aligned with Ksyp. Namely, k, could deviate from Ksyp by angle n 
(Supplementary Fig. 4). Thus, we present results for 7 < 10°, 7 < 15°, 
n< 20°, 7 < 25° and 7 < 30°. Given the marginal impact of different 
choices of n (Supplementary Figs. 5 and 6), the spectral results are 
displayed by taking the dataset for 7 < 30° as an example without loss 
of generality. The fluctuations, inthis event, count for 42% of the total 
Alfvénic fluctuations. 

To ensure the reliability of the wavenumber determination, we 
established minimum thresholds k > 1/(100d,,) and k, > 10° km“. Con- 
sequently, our observations exclude the ideal two-dimensional (2D) 
case (k, = 0). Here, T4 is infinity, indicating persistent strong nonlinear- 
ity”, Nevertheless, quasi-2D (small k,) modes are present, as k, is much 
smallerthank, at small wavenumbers (Fig. 2). These quasi-2D fluctua- 
tions satisfy T; < Ta, as shown in Fig. 3c, and exhibit weak nonlinearity. 
This weak turbulent state occurs since OB*(k,, kų)/B? is very low, 
where the Alfvénic magnetic energy density at k, and k is calculated 


with 6B (k,, kj) = as Ei e Jo Pe, (kis kyfsc) d fsc- 
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Fig. 2| Comparison between wavenumber distributions of Alfvénic magnetic 
energy Dp, (kı, k) and the theoretical energy spectra ìk, kj).a,2D 
spectral image of Dp, (k,, kj) with a high resolution (400 x 400 bins). b, 2D filled 
contours of Dp, (kı, kj) with low-resolution binning (150 x 150 bins), to clarify 
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the contours. In each panel, /,(k,, k\) at Lo = 4.6 x 10* kmis displayed by colour 
contours with black dashed curves with the same colour map as Ds, (ki, kı). The 
black dotted curves mark k = , [ki + K = 0.01/d; and 0.03/d,. These figures 
utilize the dataset for ņ < 30°. 


Evidence for the Alfvénic weak-to-strong transition 
The 2D wavenumber distributions of magnetic energy are calculated by 


De, (kuk) = f Pa, (kis kp fic) d fe. a) 
(0) 


Dg, (k,, ki) = Dg, (k1,ky)/Dg, max is normalized by the maximum mag- 
netic energy in all (k,, k,) bins, as displayed by the spectral image and 
contours in Fig. 2. Compared to the isotropic dotted curves, De, (ki, ky) 
is prominently distributed along the k, direction, suggesting a faster 
perpendicular cascade. This anisotropic behaviour is more pronounced 
at higher wavenumbers, consistent with previous simulations and 
observations”*?”?8, 

Moreover, Ds, (k,, kis compared with the modelled 2D theoreti- 
cal energy spectra based on strong turbulence*”’: 


D3 k 
o Kil ) (2) 


4/3 2/3 
Marks 


In (k,, ky) x ie exp (- 


where the injection scale Lo ~ [4.6 x 104, 8.1 x 10*] km is approximately 
estimated by the correlation time 7, = [1,300,2,300]s and 
root mean square of the perpendicular fluctuating velocity 
6V,ı = Ma cur Va = 35 km s™. InFig. 2, (ki, kı)is normalized as /,(k,, kų) 
by a constant value (one-third of the maximum magnetic energy in 
all (k,, k,) bins), as displayed by the colour contours with black 
dashed curves. The 2D distribution Dp, (k, kı) Shows two different 
properties: (1) For k; < 2 10 km, Dg, (k,,k,) is mainly concen- 
trated at k, < 7 x 10° km“ and cascading along the k, direction, sug- 
gesting that little energy cascades parallel to the background 
magnetic field, consistent with energy distributions in weak MHD 
turbulence‘. (2) For k, > 2 x 10* km", Dg, (k,, kj) starts to distribute 
to higher k,, and both wavenumber distributions and intensity 
changes of Ds, (k1, kare almost consistent with /,(k,, kų). This indi- 
cates that Dg, (ki, ky) captures some of the theoretical characteris- 
tics of strong MHD turbulence’. Besides, Dg, (ki, ky) is in good 
agreement with the Goldreich-Sridhar scaling ky « k?” (ref. 4). 
This result further confirms that the properties of Dg, (ki sky) at 
k; >2x10“*km ‘are closer to those in strong MHD turbulence. The 
change in Dz, (ki, ky) from purely stretching along the k, direction 
to following the Goldreich-Sridhar scaling k « k?’ reveals a pos- 
sible transition in the energy cascade. 


k, (km) 


Fig. 3 | Perpendicular wavenumber dependence of the compensated spectra, 
parallel wavenumber and nonlinearity parameter. a, k°/°E, N (k,) versus k,. 
The dashed line represents the scaling KP Es, (k,) « k°-?, To facilitate 
comparison with proton velocity fluctuations, magnetic field fluctuations 

are in Alfvén speed units. b, Variation of k versus k,. The dashed line represents 
the scaling kj, œ k?/°. The horizontal dotted line marks k, = 7 x 10% km”. 

c, Xp, (k., ky) spectrum calculated using the dataset for n < 30°. There are 

four zones: (1) 5 x 10° < k, < 1.6 x 10 km7, (2) 1.6 x10% < kı <3 x 10* km", 

(3)3 x10% < k, <7 x10% km^and (4)7 x 10* < k, <1x 10° km”. The first, second 
and third vertical dotted lines are around the maximum of kP’ Es, (k1), the 
beginning and the end of flattened k?’ Es, (k1), respectively. 


Figure 3a shows the compensated spectra (KŽP Ep, (k, )). The mag- 


netic energy spectral density is defined as £g, (k,) = 6B? (k,)/2ky where 
6B? (k,) is the magnetic energy density at k, (Methods). In zone 2, 
k?’ Ep, (k,)is roughly consistent with k >’? (the dashed line), indicat- 
ing that the spectral slopes of Eg, (k,) are around -2. In zone 3, on the 
other hand, kẹ’ Es, (k,)is almost flat, suggesting that Ep, (k, )satisfies 
the Kolmogorov scaling (Ep, (k1) « eo The sharp change in the 
spectral slopes of Eg, (k1) from -2 to -5/3 is apparent evidence for the 
transition of the turbulence regimes'*”. In addition, Eg, (k,) œ Kk; 
appears ina substantial portion of zone 1, indicating the weak turbu- 
lence forcing in action’. 

Figure 3b shows the variation of k, versus k, given the same Alfvé- 
nic magnetic energy. In zone 1, as k, increases, kų is relatively stable at 
k,=7x10° km. Inzone3, the variation of k, versus k agrees with the 
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Fig. 4 | The k, versus f, distributions of Alfvénic magnetic energy in the 
plasma flow frame. Ds, (ki frest) = Day (ki Srest)/DB, max is normalized by the 
maximum magnetic energy in all (k,, fies.) bins. The horizontal dotted lines 
represent theoretical Alfvén frequencies f, = |k, Val/(211), where k, = 7 x 10° km™ in 
zone 1of Fig. 3b, and k, ~ 1x 10* km in zone 1 of Supplementary Fig. 7b. The fy 
uncertainties are estimated by the standard deviation of V, (106 +11 km”), 
illustrated by error bars on corresponding horizontal dotted lines. This figure 
utilizes the dataset for 7 < 30°. 


Goldreich-Sridhar scaling k œ Ke (the dashed line). Figure 3c shows 
k, versus k, distributions of the nonlinearity parameter xg, (k1, kq), 
whichis one of the most critical parameters in distinguishing between 
weak and strong MHD turbulence’. Here, xp, (K1, ky) is calculated 
eam (Methods). For the corresponding parallel and 
ll 

perpendicular wavenumbers in Fig. 3b, Xg, (k,, kų) is much less than 
unity at most wavenumbers in zone 1, whereas in zone 3, Xp, (kK, , ky) 
increases towards unity and follows the scaling ky œ sae These results 
suggest that there is a transition from weak to strong nonlinear 
interactions, in agreement with theoretical expectations and 
simulations’*>. 

We observe a similar Alfvénic weak-to-strong transition in the 
measured proton velocity fluctuations (Supplementary Fig. 7). The 
transition scale (A,,) is estimated from the smallest perpendicular wave- 
number of strong turbulence (k, ce), where Ac = 1/k, cg. For both mag- 
netic field and proton velocity fluctuations, K; cgis around 3 x10“ km7, 
marked by the second vertical lines in Fig. 3 and Supplementary Fig. 7. 
The consistency in the transition scales estimated by magnetic field 
and proton velocity measurements further confirms the reliability of 
our findings. 

A notable perturbation is present in zone 2, as a result of local 
enhancements of the magnetic energy at k, ~1.8 x 10% km” (Fig. 2), 
leading to the simultaneous existence of strong nonlinearity (Xg, ~ 1) 
and weak nonlinearity (Xg, <1_) in the wave number range corre- 
sponding to those in Fig. 3b. Thus, the Alfvénic weak-to-strong transi- 
tion more probably occurs within a ‘region’ rather than at a critical 
wavenumber. We do not discuss the fluctuations in zone 4. The devia- 
tions in the data for 7 < 10° and ņ < 15° in zone 3 of Fig. 3b are probably 
dueto the limited number of data samples (Supplementary Fig. 6). The 
uncertainties mentioned above do not affect our main conclusions. 

Figure 4 presents f.e versus k, distributions of the magnetic 
energy, where fex is the frequency in the plasma flow frame. At 
k, <5x10° km‘, the magnetic energy is concentrated at fes: = fa, where 

Ff, isthe Alfvén frequency (horizontal dotted lines with error bars). At 
k,>1x10* km‘, the range of fex broadens, mostly deviating from f}. 
Nevertheless, the boundary of fluctuating frequencies is roughly con- 
sistent with the scaling frest « k?’ (the dashed line), indicating that 
the magnetic energy at these wavenumbers satisfies the scaling 
ky œ k?’ due to fes% ky for Alfvén modes. These results suggest that 
Alfvénic fluctuations with strong nonlinear interactions do not agree 


from 


Table 1 | Transition wavenumbers are determined by 
magnetic field measurements 


Weak MHD turbulence Strong MHD turbulence 
k versus k, Purely perpendicular Goldreich-Sridhar cascade, 
distributions cascade, k,<2x10 km”  k,>2x10 km" 
of magnetic 
energy 


Spectralslopes Wave-like (-2), Kolmogorov-like (-5/3), 


of magnetic 1.6x10%<k,<3x10%km' 3x10<k, <7x10 km" 
energy 

Nonlinearity Xea <1, k, <1x10 km" Xaa=1 and Xpq>1, 

parameter (Yq) 3x104<k,<7x10%km" 
Frequency- Single-frequency Broadening-frequency 
wavenumber fluctuations, f,est fy fluctuations with frest % Kee 


distributions k,<5x10°km"™ 


boundary, k,>1x10° km” 


with linear dispersion relations but satisfy the wavenumber scaling of 
Alfvén modes. The change from single-frequency to broadening- 
frequency fluctuations with increasing k, suggests that there is a pos- 
sible transition of turbulence regimes. 


Discussion 

The Alfvénic transition of weak to strong turbulence during cascades 
to smaller scales is one of the cornerstones of the modern MHD theory. 
Despite being proposed decades ago, evidence that confirms the exist- 
ence of the Alfvénic transition is lacking. In this paper, we present direct 
evidence of the Alfvénic transition at different angles, for example, 
the transition of energy spectra (Fig. 3a), the Goldreich-Sridhar-type 
envelope for the nonlinear parameter (Fig. 3c) and the spread of fest 
on small scales (Fig. 4; see Table 1 for a summary). Our observation 
demonstrates that the Alfvénic transition to strong turbulence is bound 
to occur with an increase of the nonlinearity, even if fluctuations on 
large scales are considered as ‘small amplitude’ (M; turb ~ 0.33). Note 
that the plasma parameters in the analysed event are generic and that 
the Alfvénic weak-to-strong transition can occur in other astrophysical 
and space plasma systems. The impact of our findings goes beyond the 
study of turbulence to particle transport and acceleration™”, magnetic 
reconnection®’™, star formation” and all other relevant fields. See, 
for example, refs. 37,38. 


Methods 

GSE coordinates 

We use the GSE coordinates in this study. X,<; points towards the Sun 
fromthe Earth, Z,,, orients along the ecliptic north pole and Y,,, com- 
pletes a right-handed system. 


Trace power spectral densities 

The trace power spectral densities of the magnetic field and proton 
velocities (Pg = Psx + Pay + Pg,z and Py = Pyx + Pyy + Pyz) were calcu- 
lated by applying a fast Fourier transform with three-point centred 
smoothing in GSE coordinates. We chose the intermediate instant of 
each time windowas the time point at which the spectral slope varies 
with time. 


Alfvén mode decomposition 

We calculated the wavenumber-frequency distributions of the Alfvénic 
magnetic field and proton velocity power with an improved Alfvén 
mode decomposition method. This method combines the linear 
decomposition method’, singular value decomposition (SVD)” and 
a multispacecraft timing analysis”. We performed the calculations 
in each moving time window with a length of 5 h and a moving step 
of 5 min. The window length selection (5 h) provides low-frequency 
(large-scale) measurements while ensuring that B, approaches the 
local background magnetic field. 
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First, we obtained the wavelet coefficients (W) of the magnetic 
field and proton velocity using Morlet-wavelet transforms”. To elimi- 
nate the edge effect due to the finite length of the time series, we per- 
formed the wavelet transforms for twice the time window length and 
cutoff of the periods affected. 

Second, the wavevector directions (K.yp(t, f.)) were determined 
by SVD for the magnetic wavelet coefficients”. SVD was used to create 
areal matrix equation (S - ksyp = 0) equivalent to the linearized Gauss’s 
law for magnetism (B - ksvp = 0). Notice that the minimum singular 
value of the real matrix S$ (6 x 3) is the best estimate of wavevector 
directions but cannot be used to determine the wavenumbers. Since 
the relative satellite separations were much shorter than the 
half-wavelength of MHD scales, the properties of the fluctuations 
simultaneously measured by the four Cluster spacecraft are similar. 
Thus, the average wavevector direction and background magnetic field 
are given by ksv = + Dic23,4 Ksvo,c,and Bo = $ Dina23,4 Bo,c, C:denotes 
the four Cluster spacecraft. 

Third, we extracted Alfvénic components from proton velocity 
fluctuations based on their incompressibility (Ksyp - 6V, = 0) and per- 
pendicular fluctuating directions (bo - 6V, = 0) in wavevector space, 
where 8V, is expressed by vectors of velocity wavelet coefficients, 
Ksyp = Ksyp/|Ksyp|, and by = Bo/|Bo|. Similarly, Alfvénic magnetic field 
fluctuations were extracted from ksyp - 6B = O and by - 6B = 0, accord- 
ing to the linearized induction equation 


W6B = k x (Bo x 6V,) Iki Ksyp x (Bo x Vp), (3) 


where kis the wavevector. Thus, Alfvénic proton velocity and magnetic 
field fluctuations are in the same direction ksyp x Bo/|Ksyp x Bo|(See the 
schematic in Supplementary Fig. 4). 

Fourth, the Alfvénic magnetic power at each time t and f,. was 
calculated as Px, (t,fec) = ; Diei23,4 Wena Wana The Alfvénic proton 
velocity power was calculated as Py, (tfs) = Wy,.c, Wi es This is 
because magnetic field data were available from the four Cluster 
spacecraft, whereas proton plasma data were available only from 
Cluster-1 during the period analysed. 

Fifth, as SVD does not give the magnitude of the wavevectors, we 
calculated wavevectors (k,(¢, f,,)) using the multispacecraft timing 
analysis based on phase differences between the Alfvénic magnetic 
field from the four spacecraft”. The magnetic field data were interpo- 
lated to a uniform time resolution of eight samples per second to give 
sufficient time resolution. We consider that the wavefront was moving 
in the direction fh with velocity V,. The wavevector k, = 27 f.m, where 
the vector m = fi/V, and the subscript A represents the Alfvénic 
component. 


f-ri Oty 
r-r; |m=|] ôt |, (4) 
r4- ri Ot, 


where Cluster-1 has arbitrarily been taken as the reference. The 
left-hand side of equation (4) contains the relative spacecraft separa- 
tions. The right-hand side of equation (4) represents the weighted 
average time delays, estimated by the ratio of six phase differences, 
j= arctan(S(Ws ), ROWE )) , to the angular frequencies, w, = 21 f.., 
where @, is from all spacecraft pairs (ij = 12, 13, 14, 23, 24, 34)). sand R 
are the imaginary and real parts of the cross-correlation coefficients, 
respectively. The four Cluster spacecraft provide six cross-correlation 
coefficients”, that is, Wye = (Wena Wg o) Wi = (Wena We. a)» Wie = 
(We, .c, Wg c?) 7 Ws = (We, c, Wac?) A We = (Wp, .c, Wac?) and W = 
(Wences Wi, c) where the angular brackets denote a time average over 
256 sto ensure the reliability of the phase differences. 


Note that a timing analysis was used to determine the actual 
wavevectors of the Alfvénic magnetic field. In contrast, SVD determines 
the best estimate of the wavevector sum in three magnetic field com- 
ponents”. Thus, k, is not completely aligned with kvp. Besides, we 
restrict our analysis to fluctuations with small angle 7 between ksypand 
k, to ensure the reliability of the extraction process (the third step). By 
relaxing the 7 constraints, more sampling points are used so that the 
uncertainty from the limited measurements decreases. On the other 
hand, with relaxed 7 constraints, k, deviates more from Kgyp, which may 
increase the uncertainty. This letter presents results from five datasets 
forn <10°,ņ < 15°, <20°,n< 25° and n < 30° to investigate the effects 
of uncertainties introduced by the combination of SVD and a 
timing analysis. 

Sixth, we constructed a set of 400 x 400 x 400 bins to obtain 

wavenumber-frequency distributions of the magnetic power 
Ps, (kı, Kp fsc) and proton velocity power Py, (k,,k).foc), where the 
parallel wavenumber is kų = ką - band the perpendicular wavenumber 
isk, =,/k} — i Each bin subtended approximately the same k,,k, and 
f To cover all MHD wavenumbers and ensure measurement reliability, 
we restricted our analysis to fluctuations with 1/(100d,.) < k < min(0.1/ 
max(d,, rai), 1/d,.) and 2/t* < frest <f,;/2, and fluctuations beyond these 
wavenumber and frequency ranges were set to zero. Here, d,, is the 
relative satellite separation, min(*) and max(*) are the minimum and 
maximum, d,is the proton inertial length, r,,isthe proton gyro-radius, 
t* is the duration studied, fest =fec — Ka - V,/(21) is the frequency in the 
plasma flow frame and V, is the proton bulk velocity with the spacecraft 
velocity being negligible. This study utilizes the representation of 
absolute frequencies: 


(frest» ka), fres > 0, 
(frest Ka) = (5) 
(frest —Ka), fes <0. 


Pe, (ki, ky .foc) were obtained by averaging P,, (k, kjfsc t) over effective 
time points in all time windows at each f.. and each k, where € = Vor B 
represents the proton velocity (V) and magnetic field (B). 


Alfvén speed units 
For comparison, this study presents the fluctuating magnetic field in 
Alfvén speed units, which are normalized by , /z9m,No, where pois the 


vacuum permeability, m, the proton mass and N, the mean proton 
density. 


Magnetic energy spectral density 
For this study, we define the energy spectral density of the magnetic 


field as Fg, (k,) = ; a where the Alfvénic magnetic energy density 
is calculated as 6B(k,) = 2 yo So Pe, (kis ky foc) d fsc- 


Nonlinearity parameter 
The nonlinearity parameter was estimated with yp, (kik) © 


k, 6By(k,, ky)/(kBo), where the Alfvénic magnetic energy density was 
calculated as 6B2 (k,, ky) = ee D =k Jo Pe, (kas ky foc) d foo Where Bo 
in Alfvén speed units was around 106 km s”. 
Frequency-wavenumber distribution of the magnetic energy 
The frequency-wavenumber distributions of the magnetic energy 
no Pay (KiKi fee)Mfec 


ky=0 
and were transformed into the plasma flow frame by correcting for the 
Doppler shift fest =f- Ka - V/(2T1). 


were approximately estimated by Dg, (kif) ¥ È 


Data availability 


The Cluster data are available at https://cdaweb.gsfc.nasa.gov. 
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Code availability 
Data analysis was performed using the IRFU-MATLAB analysis package 
available at https://github.com/irfu/irfu-matlab. 
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